Scanning Tunneling Microscopy /Spectroscopy of Vortices in LiFeAs 
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We investigate vortices in LiFeAs using scanning tunneling microscopy /spectroscopy. Zero- field 
tunneling spectra show two superconducting gaps without detectable spectral weight near the Fermi 
energy, evidencing fully-gapped multi-band superconductivity. We image vortices in a wide field 
range from 0.1 T to 11 T by mapping the tunneling conductance at the Fermi energy. A quasi- 
hexagonal vortex lattice at low field contains domain boundaries which consist of alternating vortices 
with unusual coordination numbers of 5 and 7. With increasing field, the domain boundaries become 
ill-defined, resulting in a uniformly disordered vortex matter. Tunneling spectra taken at the vortex 
center are characterized by a sharp peak just below the Fermi energy, apparently violating particle- 
hole symmetry. The image of each vortex shows energy-dependent 4-fold anisotropy which may 
be associated with the anisotropy of the Fermi surface. The vortex radius shrinks with decreasing 
temperature and becomes smaller than the coherence length estimated from the upper critical field. 
This is direct evidence of the Kramer-Pesch effect expected in a clean superconductor. 

PACS numbers: 74.70.Xa, 74.55,+v, 74.25.Uv, 74.20.Rp 



I. INTRODUCTION 

Vortices in type-II superconductors have attracted 
much attention. Since superconductivity is destroyed 
at the vortex center, the vortex core can be regarded 
as a potential well surrounded by the superconducting 
gap (SG). Therefore, quasi-particle states in the vortex 
core provide us with useful information of the SG. In 
this respect, it is particularly interesting to study vor- 
tices in iron-based superconductors [l|-Q. Characteris- 
tic features of iron-based superconductors, namely high 
transition temperature, multi-band nature, [5[ small car- 
rier density and possible sign-reversal in the SG, 
should bring about unusual quasi-particle states in and 
around the vortex core. Here we investigate vortices of 
an iron-based superconductor using scanning tunneling 
microscopy / spectroscopy (STM / STS) . 

In general, the excitation spectrum in the vortex core 
consists of discrete bound states separated by ~ A 2 /Ep, 
where A is the SG and Ep is the Fermi energy. [15j In 
conventional superconductors, Ep is many orders of mag- 
nitude larger than A and the excitation spectrum can be 
regarded as continuous. Depending on the quasi-particle 
lifetime r associated with impurity scattering, the vortex 
core can be classified into either dirty or clean regimes. 
If the lifetime broadening h/r is greater than A (equiv- 
alent to £ > Z, where £ is the coherence length and I is 
the mean- free path), quasi-particles are frequently scat- 
tered in the core and the core is well approximated as a 
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normal metal. The influence of SG on the core states is 
masked in this dirty case. On the other hand, if h/r < A 
or £ < i, quasi-particles can travel from one side of the 
vortex core to the other without being scattered. The 
most prominent feature of such a clean vortex core is 
that local density-of-states (LDOS) at the vortex center 
exhibits a peak at Ep because low-energy bound states 
are localized near the center. This LDOS peak was first 
observed by Hess et al in NbSe2 [16| and subsequent 
theoretical analysis revealed that the energy and spatial 
evolution of the LDOS peak reflect the nature of SG as 
well as normal-state band structure [l7l-[2Cij. 

If the ratio Ep/A, which roughly corresponds to the 
number of bound states, becomes small, energy separa- 
tion between the bound states is large and the excita- 
tion spectrum is no longer continuous, provided h/r < 
A 2 /Ep. Unfortunately, this so-called quantum-limit 
regime may not be reachable in conventional phonon- 
mediated superconductors. However, in iron-based su- 
perconductors, because of small carrier density (which 
leads small Ep) and high transition temperature (which 
means large A), the quantum-limit regime may be an- 
ticipated. Theoretical analysis of the quantum-limit vor- 
tex predicts strong particle-hole asymmetry in an LDOS 
spectrum at the vortex center and a Friedel-like oscilla- 
tion in the spatial variation of LDOS. [2l| 

So far, vortices are successfully observed by STM/STS 
in some iron-based supercondcutors. [22l-[24j The LDOS 
peak associated with the bound states has been detected 
in the tunneling spectra at the vortex centers of FeSe 
and Bao.6Ko.4Fe2As2. [23|, [24| In the latter compound, 
possible quantum-limit behavior is discussed based on 
the asymmetry of the spectrum. [23| However, detailed 
electronic structure of the vortex core and effects of multi 
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band are still elusive. 

In addition to the quasi-particle states in the vortex 
core, properties of many vortices, vortex matter, have 
provided a fertile area of research. In general, vor- 
tices tend to form a triangular lattice (Abrikosov lattice) 
due to the repulsive interaction between them. How- 
ever, the Abrikosov lattice may qualitatively be altered 
by quenched disorders and/or thermal agitation. [25| 
The competition among associated energy scales, namely, 
elastic energy of the Abrikosov lattice, pinning energies 
of quenched disorders and thermal energy, brings about 
various vortex phases: hexatic Abrikosov- like phase (elas- 
tic energy dominant), vortex glass or amorphous vortices 
(pinning energy dominant) and vortex liquid (thermal en- 
ergy dominant). The effect of thermal agitation and con- 
sequent vortex-lattice melting transition has extensively 
been studied in cuprate high-temperature superconduc- 
tors. However, competition between elastic and pinning 
energies is still to be clarified. 

The stability of the Abrikosov lattice is determined 
by its shear modulus Cqq which is roughly given by 
B 2 2 b(l — b) 2 where B C 2 is the upper critical field and b is 
the magnetic field B normalized by B C 2. [26| At low and 
high b where Cqq is small, the Abrikosov lattice should 
be destroyed by the quenched disorders and the amor- 
phous vortex matter may emerge. The vortex lattice 
partially retains its order in the intermediate b range but 
should contain dislocations. [27j Numerical simulation of 
the vortex configuration shows that vortices with unusual 
coordination numbers of 5 and 7 are aligned alternately, 
forming domain boundaries. [28| A crossover to the high- 
field amorphous phase occurs through roughening of the 
boundaries. 

In this article, we report results of comprehensive 
STM/STS experiments on vortices in an iron-based su- 
perconductor LiFeAs. LiFeAs is one of the very few sto- 
ichiometric iron-based superconductors and thus clean 
single crystals are available. In addition to this, cleaved 
surfaces studied by STM/STS are electronically neutral, 
which is highly advantageous to minimize the surface- 
related extrinsic effects. Tunneling spectrum exhibits 
clear two-gap features indicating that the SG amplitude 
is different depending on the Fermi-surface cylinder. In 
contrast to other vortex-imaging techniques such as Bit- 
ter decoration, STM/STS can image vortices even at high 
magnetic field where B is almost uniform. By taking this 
advantage, we study the structure of vortex matter over 
wide B range and find a signature of high-field crossover 
to the amorphous phase. We also examined the electronic 
state of the vortex cores. Tunneling spectra at the vor- 
tex center are strongly particle-hole asymmetric with a 
sharp LDOS peak just below Ep- We compare the spec- 
tra and their spatial evolution with those expected in the 
quantum-limit vortex core and show that the LDOS fea- 
tures associated with the larger SG are consistent with 
the quantum-limit behavior, but those for the smaller SG 
including the sharp peak just below Ep are not compat- 
ible with the quantum-limit vortex-core state. We also 



found that vortex radius shrinks with decreasing tem- 
perature. This is a direct evidence of so-called Kramer- 
Pesch effect [29|,[30| associated with the low-energy bound 
states in the vortex core. 



II. EXPERIMENTAL METHOD 

LiFeAs single crystals were grown by a LiAs self-flux 
technique. [31( A resistivity p ratio p(300 K) / p(20 K) of 
as high as 45 and no detectable residual specific heat 
coefficient are observed. These data demonstrate the 
high quality of the crystals. Superconducting transition 
temperature T c defined at zero-resistance is about 18 K. 
Since LiFeAs degrades in an ambient atmosphere, all of 
the sample handling procedures were carried out inside 
of glove boxes filled with purified Ar or N2 gases. All 
STM/STS experiments were performed in a constant- 
current mode using a commercial ultra-high- vacuum low- 
temperature STM system modified by ourselves (Unisoku 
USM-1300 with Nanonis controller). [32| Tunneling con- 
ductance, which is proportional to the LDOS at given 
bias voltage, was measured by standard lock-in technique 
using the built-in software lock-in detector of Nanonis 
controller. Clean and flat surfaces perpendicular to the 
[001] axis were prepared by in situ cleaving at liquid N2 
temperature. Immediately after the cleaving, samples 
were transferred to the STM unit which was precooled 
below 10 K. Electrochemically etched W wires, which 
were cleaned, sharpened and characterized by a field-ion 
microscope, were used for scanning tips. Magnetic fields 
were applied along the [001] axis. Whenever magnetic 
field was changed, the sample was heated up above T c in 
advance and then it was field-cooled to the measurement 
temperatures in order to ensure the homogeneous field 
distribution inside the sample. 



III. RESULTS AND DISCUSSION 

A. Topography and spectra in zero magnetic field 

We show a typical STM topography in Fig. 1(a). Al- 
though various types of natural defects exist, most of the 
cleaved surface is clean and flat enough to image a reg- 
ular square atomic lattice as shown in the inset. The 
lattice constant is determined to be 0.39 nm, indicating 
that the imaged atomic site is either the As site or Li 
site. 

Figure 1(b) represents the tunneling spectrum mea- 
sured in the clean region. The spectrum is skewed V- 
shaped bottomed at Ep with a bump at about 40 meV. 
We confirmed that the spectrum does not exhibit obvi- 
ous spatial variations except in the vicinity of defects. 
Results of spectroscopic studies near defects will be pub- 
lished everywhere. 

At low energy, SG structure is clearly identified. As 
shown in Fig. 1(c), spectral weight over a finite energy 
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FIG. 1: (color online), (a) STM topograph of cleaved (001) surface of LiFeAs at 0.54 K. The set-up conditions for imaging were 
sample-bias voltage V s = +20 mV and tunneling current It = 10 pA. The inset shows a magnified topography with atomic 
resolution (V s = +20 mV, It = 100 pA). (b) Tunneling spectrum taken away from defects. Tip stabilization condition was 
V s = +98.9 mV and It = 100 pA. Measurement was done with a bias modulation amplitude V mo d = lmVrms- Data were taken 
at 1.5 K. (c) Low-energy tunneling spectra at different temperatures. Two gaps are identified in the spectra and both gaps 
disappear above 18 K which corresponds to bulk T c . Tip was stabilized at V s = +25 mV and It = 125 pA. Vmod = 0.1 mVrms- 



range near Ep is completely removed at low tempera- 
tures. This is a direct evidence that the SG is finite 
everywhere on the Fermi surface. Two pairs of peaks 
can be clearly identified, implying that the SG consists 
of at least two components. By defining the SG am- 
plitude 2A using the peak-to-peak separation, the SG 
ratio 2A//c#T c (ks denotes Boltzmann constant) for the 
smaller SG (A <~ 2.9 meV) is estimated to be 3.7, which 
is close to the value expected from the weak-coupling 
Bardeen- Cooper- Schrieffer theory. The SG ratio for the 
larger SG (A ~ 6.0 meV) takes a very large value of 
7.7. In both SGs, the rising edges at lower energy 
have finite width which is wider than our energy reso- 
lution (~ 0.1 meV). This result is consistent with the 
anisotropic SG suggested by recent angle-resolved pho- 
toemission spectroscopy (ARPES) experiments [33l435j 
and quasi-particle interference effect. [36| 

The observed SG amplitudes are quantitatively consis- 
tent with the ARPES data. (33^351 We identify that the 
smaller SG opens on the outer-most hole cylinder and the 
larger SG is associated with the inner hole cylinders. A 
contribution from the electron cylinders, where ARPES 
observed A ~ 4 meV, [33|, [34[ is not clear in the tunnel- 
ing spectrum. Possible reason for this apparent absence 
may be that the gap-edge peak for the electron cylinders 
overlaps with those for the hole cylinders, or electron 
cylinders are "masked" by the tunneling matrix element. 

At elevated temperatures, both gaps are smeared out 
and disappear at about bulk T c of 18 K, indicating that 
they faithfully represent bulk superconductivity. 



B. Imaging vortex matter 

Since superconductivity is destroyed at the vortex cen- 
ter, LDOS below the SG energy should be enhanced 
around the vortices. Therefore, vortices can be imaged 



by spectroscopic imaging technique in which tunneling 
conductance is measured at every pixel of STM image. 
Figures 2(a)-(c) are conductance maps at Ep at rep- 
resentative £?'s, showing vortices as bright spots. We 
repeated the same experiment over wide B range from 
0.1 T to 11 T. As shown in Fig. 3(a), B dependence of 
vortex density perfectly follows the trend expected in the 
case that single vortex carries the single flux quantum of 
2.07xl0- 15 Wb. 

Interestingly, vortices do not form into the ordered 
triangular lattice. Two-dimensional auto-correlation 
(Figs. 2(g)-(i)) and Fourier-transform (Figs. 2(j)-(l)) im- 
ages are characterized by Independent ring-like struc- 
tures around the origin. At 0.5 T, weak 6-fold spots are 
identified along the ring in the auto-correlation image 
(Fig. 2(g)), indicating that hexatic correlation is devel- 
oping. At higher £?, the ring-like structures are isotropic 
in magnitude but the rings in the Fourier-transform im- 
ages are deformed into rectangular shape elongated along 
the As- As direction. The emergence of 4-fold anisotropy 
suggests that the symmetry of underlying crystal lattice 
and/or anisotropy of the SG affect the configuration of 
vortices but they are not strong enough to arrange vor- 
tices into regular square lattice which was observed in 
YN12B2C, for example. [37], [38|. Nevertheless, the ob- 
served ring-like structures indicate that there is a well- 
developed distance correlation between vortices but the 
orientational correlation is almost absent. 

Such a highly disordered vortex matter in LiFeAs is 
in sharp contrast to the vortices in Bao.6Ko.4Fe2As2 
where well-ordered triangular vortex lattice has been ob- 
served. [23| The difference between these two materials 
can be ascribed to the difference in the shear modu- 
lus C*66 associated with the difference in B C 2. [26| In 
LiFeAs, B c2 along the [001] axis is about 15 T, [39| 
while that of (Ba,K)Fe2As2 is estimated to be more 
than 60 T. J4(J, [4JJ This difference gives rise to an 
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FIG. 2: (color online), (a)-(c) Images of vortices at 1.5 K obtained by mapping tunneling conductance at Ef- Tip was 
stabilized at V a = +20 mV and It = 100 pA. V mo d = 0.7 mV rm s. Arrows denote the nearest Fe-Fe direction, (d)-(f) Delaunay 
triangulation diagrams obtained from vortex images shown in (a)-(c). Vertices with symbols denote vortices with coordination 
numbers different from conventional value of 6. (Blue triangle: 5, red square: 7, green star: others) Dashed black square is 
drawn separately from the edge of the field of view (black solid square) by the mean value of inter- vortex distance at each field. 
Defect fraction (Fig. 3(b)) was calculated using vortices within this dashed black square to avoid the effect of the edge, (g)-(i) 
Two-dimensional autocorrelation-function calculated from vortex images shown in (a)-(c). (j)-(l) Two-dimensional Fourier 
transformations from vortex images shown in (a)-(c). 



order-of-magnitude smaller Cqq in LiFeAs than that in 
(Ba,K)Fe2As2- Therefore, vortex lattice in LiFeAs is 
much more susceptible to quenched disorders than that 
in (Ba,K)Fe2As2- In other words, LiFeAs can be a play- 
ground to investigate how the structure of vortex matter 
subject to quenched disorder changes as a function of B. 

In order to analyze the local structure of vortex matter, 
we performed Delaunay triangulation analysis as summa- 
rized in Figs. 2(d)-(f). Each vertex denotes the position 
of vortex and vertices with symbols correspond to defect 
vortices whose coordination numbers are different from 
6. We plot defect- vortex density as a function of B in 
Fig. 3(b). In the low field region, defect density increases 
rapidly and reaches more than 40 %. 

At B = 0.5 T, a domain boundary, which consists of al- 
ternating defect vortices with coordination numbers of 5 
and 7, has been observed. (Fig. 2(d)) More grain bound- 



aries are observed at 2 T but each grain boundary still 
keeps its alternating structure. (Fig. 2(e)) At 8 T, de- 
fects tend to cluster and defect vortices with extraordi- 
nary coordination numbers of 4 and 8 emerge. (Fig. 2(f)) 
In order to examine the defect structure within the clus- 
ter at 8T, we compare local vortex configurations at 2 T 
and 8 T in the fields of view which contain equal num- 
ber of vortices. (Figs. 3(c),(d)) In contrast to the domain 
boundary at 2 T, well-defined alternating line-like struc- 
ture is no longer observed at 8T. Namely, even though 
defect density saturates above ~2 T, local vortex config- 
uration still changes with B. 

These observations suggest the following picture of 
order-disorder crossover in the vortex matter. In the 
B region where the elastic energy dominates, number of 
defect vortices is small. Even if there is a defect, its co- 
ordination number should be 5 or 7, which is close to 
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FIG. 3: (color online), (a) Applied-field dependence of vortex 
density in LiFeAs. Solid black line denotes the behavior ex- 
pected in the case that single vortex carries the single vortex 
quantum of 2.07xl0~ 15 Tm 2 . (b) Applied-field dependence 
of fraction of defect vortices with coordination numbers dif- 
ferent from conventional value of 6. (c) and (d) Delaunay 
triangulation diagrams at 2 T and 8 T in the fields of view 
which contain equal number of vortices. Symbols have the 
same meanings as in Figs. 2(d)- (f). 



the regular value of 6. Such defect vortices are aligned 
alternately forming a domain boundary. This situation 
does not change so much as long as the elastic energy 
is relevant. At high enough £?, however, Cqq diminishes. 
As a result, the domain boundary becomes rough and 
the vortex matter approaches to the uniformly disordered 
amorphous phase which contains vortices with extraor- 
dinary coordination numbers less than 5 and more than 
7. These behaviors are totally consistent with the result 
of numerical simulation. 1281 



C. Shape of a single vortex 

We now focus on spectroscopic features of a single vor- 
tex. All of the data shown below are taken at B = 0.5 T 
and the reproducibility of the results has been checked 
in a different tip-sample combination. Figures 4(a)- (d) 
show conductance maps at different energies around a 
vortex. Vortex-core state at Ep exhibits a 4- fold star 
shape with high LDOS tails along the nearest As-As 
directions. Each tail splits in parallel with increasing 
energy. These behaviors closely resemble the energy- 
dependent vortex-core shape of NbSe2, [Hj] except that 
in NbSe2 vortex core has a 6-fold anisotropy rather than 
4-fold. 



Theoretical models exist which account for the energy- 
dependent vortex-core shape include effects of vortex lat- 
tice [13, G3, SG anisotropy pjl [20| and crystal-lattice 
potential which determines the shape of the Fermi sur- 
face. [nlHIIl 

In the case of LiFeAs, the effect of vortex 
lattice can be safely excluded because vortices at 0.5 T 
form into a distorted hexatic lattice (Fig. 2(a)), which 
does not fit with the observed 4-fold anisotropy of the 
vortex core. 

At low enough energies, LDOS distribution around a 
vortex is primarily governed by the smallest SG and the 
Fermi surface where it opens. In LiFeAs, ARPES ob- 
served the smallest SG in the outer- most hole cylinder 
and the SG exhibits 4-fold anisotropy. [33M35| The SG is 
smaller along the Fe-Fe direction than along the As-As 
direction, which would give rise to the vortex-core shape 
with the LDOS tails along the Fe-Fe direction. This is in 
contradiction to the observed LDOS distribution at Ep 
shown in Fig. 4(a), indicating that the SG anisotropy 
does not play a major role. 

On the other hand, the shape of the Fermi sur- 
face is consistent with the observed anisotropy of the 
vortex-core state. According to ARPES, the outer- most 
hole cylinder has a rounded square cross section with 
flat region perpendicular to the As-As direction. (43| 
This anisotropy should result in a larger number of 
quasi-particles traveling along the As-As direction, form- 
ing high LDOS tails around a vortex. Indeed, Wang, 
Hirschfeld and Vekhter analyzed the LDOS distribution 
around a single vortex in LiFeAs with the assumptions of 
various SG anisotropy and concluded that Fermi-surface 
anisotropy in LiFeAs is strong enough to dominate the 
various types of SG anisotropics. [44] 



D. LDOS spectrum and its spatial evolution 

In order to get more insight into the vortex core state, 
we study the individual tunneling spectrum and its de- 
tailed spatial evolution. Figure 4(e) shows the spectrum 
taken at the vortex center. We also plot the spectrum 
away from vortices where the two-gap feature observed 
in the absence of B is well maintained. The spectrum at 
the center exhibits a pronounced peak just below Ep (- 
0.9 meV) as well as some other weaker peaks. False color 
plots of spatial evolutions of tunneling spectrum along 
Fe-Fe and As-As directions are shown in Fig. 4(f) and 
(g), respectively. The LDOS peak at ~-0.9 mV shifts 
further away from Ep with increasing distance from the 
vortex center and smoothly approaches the smaller SG. 
The counter branch is also observed at symmetric en- 
ergies above Ep, although the intensity is much weaker 
near the vortex center. There is an energy-gap-like fea- 
ture apparent at Ep even at the vortex center. Spectral 
weight at low energies is higher along the As-As direc- 
tion than along the Fe-Fe direction, giving rise to the 
anisotropic shape of the vortex. 

LDOS features associated with the larger SG show 
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FIG. 4: (color online), (a)-(d) Tunneling-conductance images 
around a single vortex at different energies. Arrows denote the 
nearest Fe-Fe direction. Tip was stabilized at V s = +25 mV 
and I t = 250 pA. V mo d = 0.2 mV rm s- Data were taken at 
0.54 K. (e) Tunneling spectra taken at the center of vortex 
(red) and away from vortices (blue). V mo d = 0.10 mV rm s. 
Line profiles of tunneling conductance across the vortex cen- 
ter, (f) Along the nearest Fe-Fe direction, (g) Along the 
nearest As-As direction. V mo d = 0.15 mV rms . 

more complicated behavior. A broad LDOS peak at ~- 
2.3 mV at the vortex center exhibits jumps in its energy 
with increasing distance from the center and finally de- 
velops as the peak at the larger- SG energy. Above Ep, 
the gap-edge peak evolves from the peak at ~+4.2 mV 
at the center. Namely, the quasi-particle branches asso- 
ciated with the larger SG are asymmetric between filled 
and empty states not only in their intensity but also in 
their energy. These features are in sharp contrast to the 
smooth particle-hole symmetric evolutions observed for 
the quasi-particle branches associated with the smaller 
SG. 



The observed vortex core states in LiFeAs are appar- 
ently different from those of conventional superconduc- 
tors, e.g., NbSe2- In NbSe2 where the vortex-core ex- 
citation spectrum in almost continuous, LDOS peak is 
located just at Ep- [16] With increasing distance from 
the vortex center, the peak splits along smooth branches 
which are symmetric in both intensity and energy be- 
tween filled and empty states. [16| By contrast, in LiFeAs, 
spectral weight in the vicinity of vortex center is strongly 
particle- hole asymmetric, there is a gap-like feature even 
at the vortex center, and evolution of quasi-particle 
branches for larger SG is neither smooth nor symmetric. 
These observations indicate that factors which are lack- 
ing in the vortex of conventional superconductor play a 
vital role. 

Proximity to the quantum limit is a plausible candi- 
date of the observed asymmetry, as has been suggested in 
Bao.6Ko.4Fe2As2. [23] In the quantum-limit vortex, the- 
oretical calculation of LDOS spectrum predicts strong 
particle-hole asymmetry at the vortex center which is as- 
sociated with the lowest bound state. [21| LDOS peak 
for the lowest bound state appears only on e mpty (filled) 
state if the carriers are electron (hole). [U, [23j This is 
apparently consistent with the observed asymmetry in 
LiFeAs; the LDOS peaks are located below Ef and the 
relevant SGs are opening in the hole cylinders. 

In this scenario, energy of the LDOS peak \E P \ at 
the vortex center represents the energy of the lowest 
bound state given by A 2 /2Ef- For the larger SG, we 
observed that E p ~ —2.3 meV and A ~ 6.0 meV, lead- 
ing Ef ~ 7.8 meV. This Ef value is very small, in gen- 
eral sense, but does not contradict with the ARPES re- 
sults (3314351. l43j where the top of the inner hole band 
is very close to Ef- The particle- hole asymmetric spa- 
tial evolution of the LDOS peak, especially the step-like 
change in its energy in the filled state, may be due to the 
Friedel-like oscillation. [2l| Namely, it is possible that 
the observed asymmetric quasi-particle states associated 
with the larger SG reflect the quantum nature of the vor- 
tex electronic states in LiFeAs. 



However, the quantum-limit picture may not apply 
for the smaller SG because calculated E F ~ 4.7 meV 
from E p ~ —0.9 meV and A ~ 2.9 meV is much 
(at least order-of-magnitude) smaller than the value ex- 
pected from ARPES dispersion. [3314351. l43j In addition, 
the quasi-particle branch approaching to the smaller SG 
is observed to be very smooth without any trace of 
Friedel-like oscillation. Such a smooth evolution may 
suggest that quantum nature is absent in the vortex elec- 
tronic states associated with the smaller SG. It is an im- 
portant open issue to identify what gives rise to the pro- 
nounced particle-hole asymmetry in the spectral weight 
and opens an energy gap at the vortex center. 
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FIG. 5: (color online). Temperature dependence of the tun- 
neling spectra at the vortex center (a) and away from vortices 
(b). Tip was stabilized at V 8 = +25 mV and I t = 250 pA. 
Vmod = 0.1 mVrms- Each curve (except at 7.5 K) is shifted 
vertically for clarity. Thin black lines denote the convoluted 
0.54 K spectra by Fermi-Dirac distribution function. 



E. Temperature dependence: Kramer-Pesch effect 

Next we examine temperature T dependence of the 
vortex-core states. Figure 5(a) shows the tunneling spec- 
tra taken at the vortex center at different T. The LDOS 
peak at -0.9 meV diminishes rapidly with increasing T. 
We compare the observed spectrum at each T with the 
spectrum at 0.54 K convoluted with the Fermi-Dirac dis- 
tribution function at T. It is clear that the observed 
LDOS peak is broader than that expected from the ther- 
mal smearing effect. This is in stark contrast to the T 
dependence of the spectrum away from vortices where 
each spectrum is well reproduced by the convoluted spec- 
trum as shown in Fig. 5(b). This contrasting behavior 
strongly suggests that the temperature not only smears 
the quasi-particle states but also changes the excitation 
spectrum in the vortex core. 

The most plausible origin of this unusual temperature 
effect is the so-called Kramer-Pesch (KP) effect [29l |30|. 
In clean superconductors where h/r < A, quasi-particle 
distribution below A depends on T. The lower-energy 
bound states are more and more populated at lower T. 
This T-dependent quasi-particle distribution gives rise 
to the T dependence of spatial variation of SG around 
a vortex, and vice versa. Kramer and Pesch first ana- 
lyzed this process using quasi-classical Eilenberger equa- 
tion and pointed out that, with lowering T, the spatial 
variation of SG deviates from ~ tanh(r/£) behavior ex- 
pected from Ginzburg-Landau theory. [29j Here, r de- 
notes distance from the vortex center. In this case, it is 
convenient to introduce a new length scale £i defined as 
£i = lim r ^o(^/A(r))A(r -+ oo) to measure the vortex- 
core size. Kramer and Pesch found that £i shrinks lin- 
ear in T, provided k B T > A 2 /E F [29|. Later on, Gygi 
and Schliiter analyzed a clean vortex core based on the 
quantum Bogoliubov-de Gennes formalism. As well as 
reproducing the shrinkage of vortex core, they clarified 
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FIG. 6: (color online). Tunneling-conductance images at Ef 
at different temperatures. Tip was stabilized at V s — +25 mV 
and I t = 250 pA. V mo d = 0.2 mV rm s. Note that color scale is 
different from that for Fig.4(a)-(d). 



the T dependence of vortex bound states. [18[ 

For the sake of examining the relevance of the KP ef- 
fect, the analysis of the T dependence of tunneling spec- 
trum would not be appropriate because individual bound 
states are not resolved in Fig. 5(a). Instead, the T de- 
pendence of the vortex-core size can be a direct evidence 
of the KP effect. Figure 6 shows the conductance maps 
at Ep taken at various T. It is clear that apparent vortex 
size increases with increasing T, as expected from the KP 
effect. In order to quantify this process, we take conduc- 
tance profiles along Fe-Fe and As- As directions at each 
T. (Figs. 7(a) and 7(b)) We fit each profile with Gaus- 
sian function to obtain the half-width-at-half-maximum 
(HWHM), which can be a good measure of the vortex 
size. As shown in Fig. 7(c), HWHM decreases rapidly 
with decreasing T and tends to saturate below about 
2 K. We note that T dependence of £ estimated from 
B C 2 [39[ is very weak in the T range studied and does 
not account for the observed strong T dependence. In 
addition, HWHM at our lowest T= 0.54 K (1.9 nm and 
2.5 nm along Fe-Fe and As- As directions, respectively) 
is much shorter than £(T = 0) ~ 4.64 nm. [39[ These 
observations strongly suggest that the KP effect plays an 
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FIG. 7: (color online). Line profiles of tunneling-conductance at Ef at different temperatures, (a) Along the nearest Fe-Fe 
direction, (b) Along the nearest As- As direction. Solid lines denote results of Gaussian fitting, (c) Temperature dependence of 
half-width-at-half-maximum of the fitted Gaussian peak, which represent the vortex-core size. 



important role in LiFeAs. 

The importance of the observation of KP effect in 
LiFeAs is twofold. Firstly, to our knowledge, the KP 
effect has only been detected by muon spin rotation tech- 
nique. [45| Muon spin rotation primarily senses the field 
or current distribution inside the sample and the vortex- 
core size is defined as the distance from the vortex cen- 
ter at which the current density reaches the maximum 
value. [Ill In contrast to this rather indirect way which 
requires a specific model, STM/STS gives a direct unam- 
biguous evidence of KP effect. 

Secondly, the observation of KP effect and its satu- 
ration at low T allows us to estimate the energy of the 
lowest bound state. In the quantum- limit vortex, vortex 
core stops shrinking at T < A 2 /ksEp^ because no more 
low-energy states are available. [211, [29j In the presence of 
impurity scattering, the saturation occurs even at higher 
T. [30| The observed saturation at T ~ 2 K sets the up- 
per bound of the lowest bound-state energy associated 
with the smaller SG to be about 0.2 meV. 



IV. CONCLUSIONS 

Comprehensive STM/STS experiments on LiFeAs over 
wide B and T ranges provide new insights into the prop- 
erties of vortex matter and the electronic states of indi- 
vidual vortices. 

Zero-field tunneling spectra contain two superconduct- 
ing gaps. By comparing with ARPES results, [33l-[35| we 
identify the larger (A ~ 6.0 meV) and smaller (A ~ 
2.9 meV) gaps as opening on inner and outer hole Fermi 
surface cylinders, respectively. 

Vortices in LiFeAs do not form into a well-ordered lat- 
tice while there is a well-developed distance correlation 
between vortices. We found that local hexatic correlation 
retains at low B but the vortex matter contains domain 
boundaries which consist of alternating vortices with un- 
usual coordination numbers of 5 and 7. With increasing 



B, the vortex matter approaches to uniformly disordered 
amorphous state through increasing number of domain 
boundaries and roughening of each boundary. 

LDOS distribution at Ep around a single vortex ex- 
hibits tails along the nearest As- As directions and each 
tail splits in a parallel fashion at elevated energies. This 
4-fold anisotropic shape of the vortex core is associated 
with the square shape of the outer hole Fermi-surface 
cylinder. 

The multi-gap character of LiFeAs manifests itself in 
the spatial evolution of LDOS spectrum around a vor- 
tex as two quasi-particle branches. Tunneling spec- 
trum at the vortex center is strongly particle-hole asym- 
metric, being reminiscent of the behavior expected in 
the quantum-limit vortex core where the discrete bound 
states inside the core are smeared out neither by impurity 
scattering nor thermal broadening. The quasi-particle 
branch associated with the larger superconducting gap 
exhibit step-like energy change as approaching to the vor- 
tex center. This is consistent with the Friedel-like oscilla- 
tion expected in the quantum-limit vortex. [21] However, 
quasi-particle states associated with the smaller super- 
conducting gap are not compatible with the quantum- 
limit behavior, although energy-gap like feature is ob- 
served even at the vortex center. Observation of Kramer- 
Pesch effect, strong T dependence in the vortex-core size, 
also suggests that energy-separation between the discrete 
bound states associated with the smaller gap is not large 
enough to bring about detectable energy gap at the vor- 
tex center. The mechanism which opens the energy gap 
at the vortex center is an interesting open question. Nev- 
ertheless, we anticipate that observed rich spectroscopic 
features exposed in the vortex core of LiFeAs serve as a 
touchstone of various theoretical models for the super- 
conducting state in iron-based materials. 

Note added.- During the completion of the manuscript, 
we became aware that Chi et al. also observed two su- 
perconducting gaps in LiFeAs by STM/STS Hfj. 
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